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SINGLE LINE IN SPACE SOLUTION TO GEOCENTRIC COORDINATES OF AN EARTH 
SURFACE POINT 


1. INTRODUCTION 


The usual method of determining the coordinates of a point 
on the earth's surface photogrammetrically is with intersection equations 
formed with the space camera coordinates and orientation matrices of 
an overlapping pair of photographs on which the object point is imaged. 
The usual method is the desirable method when the intersection angles 
or, sometimes referred to as parallax angles, are 10° or more. A single 
line solution requiring the space camera coordinates and the geocentric 
orientation matrix of the camera x,y,z2 axis provides a direct geometrically 
accurate solution when intersection lines enclose less than 10° angles. 
The single line solution is independent of camera cone angle geometry 


or the base-height ratios. 


Space born panoramic photography has the advantage of high 
resolution paraxial imagery and wide cone angle normal to the direction 
of motion but a limitation in an extremely narrow cone angle in the 
direction of motion. This results in weak intersection angles along the 


orbital path. 
2 GEOMETRY OF THE SINGLE LINE DETERMINATION 


Assume an orbital point o in space whose coordinates have 
been determined with time T and known orbital parameters, Also assume 
that the direction cosines of a line to an unknown point on the earth's 
surface are known by some preprocessing means to be expanded later. This 
arrangement results in a line in space piercing an ellipsoid of revolution 
from a known point o in space. The intersection of the line in space 
with the ellipsoid of revolution defines the geocentric position of 
unknown point. The geometry of the solution is illustrated in Figure l. 


The ellipsoid of revolution is expressed with the equation. 
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phe 


or 


sho 


Figure 1 GEOMETRY OF SOLUTION 
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x? + y? + Z* =1 


Where X, Y, Z are the geocentric coordinates of a surface point to be 
determined; a and b are the semi-major and minor axis of the earth. The 


length of the line in space, M, may be expressed: 


(X, - X) cos a + Y - Y) cos 8 + (25 -~ Z) cos y = M 


or 
2 2 Bee eo 
(X, X)° + (Yo y)* + (2, Z) =M 
also 
X - X =M cos a 
Y - Y=M cos B 
Zo - Z=M cos y 
or 
X =X --M cos a 
re) 
Y = Y -M cos 
re) 


Z=Z -M cos Y 
Changing the form of the ellipsoid of revolution 


(x2 + v7). p2 + 22 a? = ah? 


b? = a? - ¢? 
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or 


(x? + y?) a? (1 - e”) + Zaz = a* (1 - e”) 
Division by a? 
(x2 + v2) (-e2) + Z? =a? (1 - e*) 
Substituting 
{X,, - M cos a)? + Z, -M cos 8)7] (1 - e7) + (2, - M cos vy)? 
= a7(1- e*) 
Expanding the squares 
[M? (cos* a+ cos? 8) (l- e*) - 2M (cos OX + cos BY) (1 - e?) 
+ M’cos? y - 2M cos YZ, + Re Seige me: a es ery 2S Grane er) 0 
Collecting coefficients of descending powers of M 


ig M2[ (cos? a + cos?B) (1 - e?) + cos*y)] 
- M 2 [(cos OX, + cos BY) (1 - e”) + cos y2Z5] 
+ [&, + ve ») (lL - e?) + dee - a2(1 - e%)] =0 


which is a quadratic in M 
Let A = [(cos* & + cos*B) (1 - e7) + cos“y)] 


B = - 2 [(cos OX, + cos BY 5) (1 - e*) + cos yZ,] 
C= LX, + ee (Ll - e7) + ZS oo grt e625] 
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Then: 
ee Vv a 
M = B+ vB 4AG 
2A 
whence 


X =X -M cos a 
Y=Y -Mecos 8 
Z=2Z -Mcos y 


The quadratic yields two values for M. The obviously correct value is 
the shorter one inasmuch as the longer line is the continuation of 

the line after it pierces the first surface at the point of emergence. 
The two values of M are illustrated in Figure 2, The solution assumes 
that the elevation of the point is zero. The solution without 
modification is directly applicable to shore line points and low lying 


coastal regions at any latitude. 


An elevated point will be displaced along M away from the 
camera station toward the earth's surface. A combination of a high 
elevation and large vertical angles enclosed between the line in space 
and the camera radius vector can cause large errors in both the length 
M and the geocentric coordinates (X, Y, Z of the point if the elevation 
h is neglected. The error is not caused by neglecting h in the reduction 
of latitude and longitude but by the displacement when the point is 


projected to the earth's surface, 


Three rigorous solutions are described that yield the correct 
geocentric coordinates for a point within elevation h. A common property 
of each of the solutions is the solution for M. When h # 0 simple 
additional computations are performed to accurately correct for the 
displacement. Historically an approximate correction was made by 


differentiating the equation: 
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ge 


Figure 2 GEOMETRY OF TWO VALUES OF M 
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where 
o= (x? + vy? +27)” 


~ 22 2 2, 45 
0 Oe aan Zo) 


Y Z 
cos i = cos a an + cos B (0 + cos Y oO 


°o °o %o 


Differentiating M with respect to Pp 


20 . do = 2M .dM - 205 cos CoaM 


putting do = h 


dM= 90.h 
M- PR cos So 


The geometric significance of this equation is shown in Figure 3. The 


ratio Q = sec C and sec t.h= dM 
M- P, cos - 


Because the equation neglects earth curvature dM is over corrected 
placing P at Pi. In any case, assuming dM to be correct 
Xp X - dM cos 4 
Yp = Y - dM cos B 
Zp = Z-—- dM cos Y 
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Figure 3 APPROXIMATE ADJUSTMENT TO M FOR h 
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The three methods described remove the errors arising 


from the above approximations. 


For completeness for the case where h = 9 the geodetic 


latitude $ and longitude A are as follows: 


tan ¢ = Z 
a 
(xX? + ¥*)%° (1 - e”) 
tan A = 


a 
X 


and the geocentric latitude >' is simply tan $' = a ns 
Ce ae a) 


also for other applications 


sin ¢ = Z le 
[(x? + ¥*) (l-e7) + 27] 
cos ¢ = (X2 + y¥*) (1 = e”) : 
[m2 + v2) (1 - e2) + 2272 
sin r} = ose 
(x? + v7) 4 
cos A * x L 


(x? +. 2)" 


There is no justification for neglecting h. There are very 
few land platforms in the world that the elevation can not be interpolated 


from available maps to several hundred feet, 


It was assumed in developing the quadratic for M that the 
direction cosines of the line cos a, cos 8, cos y-and the geocentric 
coordinates of o, Xo Yo? Zo were known. The -ddtermination of’ these 
values is considered preliminary data processing. 
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The next section covers the various forms in which the orientation 
and the spatial camera coordinates may have and the procedure for 
transforming the various given forms to the necessary geocentric 


reference orientation of the camera axis and coordinates of the camera 


station. 
3, REDUCTION ‘OF GEOGENTRIC CAMERA ORIENTATION AND COORDINATES 
3.1 Geocentric Goordinates of the Camera Station 


Geocentric Loordi na 


The advantage of geocentric coordinates is the freedom 
from problems of earth curvature and local tangent planes changing 


rapidly with each position of the camera in orbit. 


The geocentric coordinates of the camera station are 
computed with Keplerian equations of an orbit. While these equations 
do not account for periodic perturbations modern orbital computing 
programs output orbital parameters at sufficiently close intervals to 
justify use of the simpler Keplerian equations. An earth centered orbiting 


satellite is illustrated in Figure 4. 


(1) Xo» Y5° Z denote an earth centered star-referenced 
acocataaes system referred to the vernal equinox in 
the plane of the earth's equator. 

(2) X 3 Yo Zi are geocentric coordinates of a point or 
camera station on the orbital path at clock time Ts 

(3)  Perigee and apogee are near and far points of the 
orbital path. These points define the long axis 
(2 a, ) of the orbit. 

(4) An elit oeieal orbit has an eccentricity eo defined 


by 
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To PERIGEE 


L- ORBITAL 


APOGEE 


Figure 4a ORBITAL PARAMETERS 
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y / 
Co | EARTH CENTER PERIGEE 
R'=Q, COS E 


Figure 4b RELATION OF CENTER OF EARTH TO CENTER OF 
ORBIT AND PARAMETRIC ANGLES E 
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cy is the distance along the long axis between the 
geometric center of the orbit and the center of the 
earth. 

(5) a, 7 is the distance from perigee to the center of 
the earth. 

(6) a + cy is the distance from apogee to the center of 
the earth. 

(7) by is the short diameter of the orbit at the geometric 
center and perpendicular to the long axis. 

(8) The path of the satellite defines a plane neglecting 
perturbations that is space oriented just as the vernal 
equinox is space oriented (fixed in space). 

(9) The intersection points of the plane of the orbit with 
the plane of the celestial equator are called nodes. 
The ascending node is the part where the satellite moves 
northward whereas the descending node is the part where 
the satellite moves southward. 

(10) 1 is the inclination of the plane of the orbit to the 
plane of the equator. 

(11) Q is the right ascension of the ascending node or the 
angle measured eastward from the vernal equinox in the 
plane of the celestial equator to the ascending node. 

(12) w is the angle subtended by the ascending node and 
perigee at the center of the earth. 

(13) v is the angle subtended by the camera at any time T 
and perigee point at time To: 

(14) ES is a parametric angle that defines the coordinates 


of a point referred to the orbital ay and De axes. 
R’ = a_cos E 
oO ° fe) 
z' = b_sinE 
fo) fe) ° 


(15) n is the mean radian velocity of an orbiting particle. 
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(16) GM is the product of a gravitational constant and the mass 
of the earth. 

(17) P-is the clock time of two successive perigee passages 
of a satellite. 


Given orbital parameters age eo» Q 4 WwW qT, and time T 


fe) ‘e) 0? 
determine Xo? Yo? and zy of the camera station referred to the vernal equinox: 


12 
where GM = 14076.92056 x 10 ft */sec*. Parametric angle Ey is computed 


first by an iterative process: 


n (T-T_) E -e sin£E 
fe) 9) oO re) 


ie t 
n (T T)) Et (approximate ED 
= aa 1 * ’ = 
n (T T) (E° e, sin ES ndT 
AE = ndt 
1 -e cos E! 
o re) 


E = E' + AE 
ce) 1e] 


Seldom are more than two iterations required. The radius 


vector D5 follows 


and then 
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2\% 
(1 - en) t/2 sin Ey 


sin Vis 
1- e cos & 
re) fe) 
2 
a l-e 
Pe oO ( 
L+e cos ) 
fe) 
also 5 
sin Eo + (1 - @5) eal 
1+e cos v 
re) 
e + cos Vv 
cos E, = 0 


1+e cos Vv 
oO 
The geocentric coordinates follow: 
with law of cosines 
cos Q = cos (vtw) cos 2 - sin (vtw) sin 2 cos 1, 
cos Bo = cos (vw) sin 2+ sin (vtw) cos 2 cos 1. 


cos y, = sin (vtw) sin t 


whence 


x, = Pg cos Oy 
Y, = P5 cos Bo 
Z, = 


P, 698 Y, 
Since the earth is rotating it may be desirable to refer the X,Y 


coordinates to the Greenwich Meridian in place of the vernal equinox. 


This is accomplished by replacing {% with 4 
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A. = GST are - 2 
fo) 
“GST = GCT + sol - sid + GST o 


GCT is the T of the exposure. In any case 


x. p, [cos (v + Ww) cos rG - sin (v + w) sin AG cos i 


— + + . 
YS plcos (v + w) sin AG sin (v + w) cos ro cos i 
Zo = Z, since Zo is the axis of rotation. 
3.2 Camera Orientation Matrix 


It is most practical to refer the camera x, y, z axes to 
the earth geocentric X, Y, Z axes. In practice the camera orientation 
is either referred to the local vertical and local meridian or the 
sidereal geocentric coordinates depending on how the camera: orientation 
is determined, If the camera orientation is determined with Greenwich 
reference geocentric object space coordinates with well-known collinearity 


equations the matrix has the form: 


x y Z 
X cos a cos cos a 
g x Z 
XY s cos cos 
g co B B B 3 
2 cos Nig cos Y. cos 5 


where 


cos &© = sini cos A +cosrA_ sin A_ sind 
x °) A p Zz Pp 


cosQ® = sind sin A_ -cosA_cos A_ sin 6 
Pp Zz Pp Z Pp 


re) 
i2) 
a 
2 

N 
i) 


cos A_ eos 6 
Pp 
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Where A, is 


meridian * 


cos A_ + sin i 
Z Pp 


cos g. = = cos i 
x Pp 
cos 8 = -cos A_ sin A: 
¥ P z 
cos 8 = sin A_ cos 6 
Z Pp 
cos y_ = sin A_ cos 6 
x Z p 
cos Y. = cos A_ cos 6 
y Zz 
cos Y_ = sin 6 
Z Pp 
the angle 


- sin Ar 
Pp 


sin A 
Zz 


cos A 
Zz 


sin 6 


sin 6 
Pp 


between the camera y axis and the apparent local 


is the apparent longitude of the optical axis translated 


parallel to itself to the center of the earth and 5 


the optical axis defines with the plane of the ae 


form of the 


orientation matrix XG 


replaces Qo 


is simply the angle 
If this is the 


in the reduction of space 


coordinate so that both orientation and coordinates refers Greenwich. 


Suppose the camera orientation is determined with star 


data exposed simultaneously with the terrain camera and having a known 


relative orientation with the star camera. 


x y 

xX cos Ox cos ay 
8 : s s 
XY cos Bx. cos BY, 
Zs cos NS cos Bh 

where cos Ox. = - sin % cos A_ + cos & 
Ss Pp Zz Pp 

cos Gy =- sin 2 sin A = cos 2 

s Pp Zz p 

17 


This is the more likely form: 


sin A 
Zz 


cos A 
Zz 


cos 


cos 


Z 


AZ 
Ss 
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cos OZ = cos as cos 6 
cos Bx = cos 2. cos A + sin. sin A_ sin 6 
s p Zz Pp Z P 
cos By. = cos 2. sin A_ - sin 2 cos A_ sin 6 
s p Zz p Zz Pp 
cos ®z = sin 2 cos 6 
Ss Pp Pp 


cos YX, = cos A, cos 6 


cos YY, sin A, cos 8 


cos YZ. = sin o 


A, and oa have the same meaning as before and ue is the 


right ascension of the optical axis. 


In still other cases the orientation matrix may be referred 
to the local vertical and the local tangent plane at the foot of the 
local vertical with object space Y lying in the local meridian. This 


form employed usually for altitudes under 10 miles has the form, 


x y Zz 
East (X) cos Ax cos ay cos OZ 
North (Y) cos Bx cos By cos Bz 
Zenith (Z) cos YX cos Yy cos YZ 


where 


cos aX = cos S cos A + sin S sin A, cos t 
cos @y = sin S cos AL - cos § sin A, cos t 
cos az = sin A. sin t 


cos 8x = -cos S sin A, + sin S cos A, cos t 
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cos By = - sin § sin A, - cos S cos A, cos t 
cos ®z = cos A, sin t 

cos yx = sin S sin t 

cos ‘Yy = cos S sin t 


cos YZ = cos t 


S, is a rotation about the optical axis from the camera y axis to the 
image of the direction of gravity. t, is the angle enclosed between the 
optical axis and the vertical and A, is the angle in the horizontal 
plane at the foot of the vertical measured clockwise from the vertical 


projection of the optical axis to the Y axis lying in the local meridian, 


The set of Eulerian angles defining each matrix is 


illustrated in Figures 5a, 5b, 


The basic point of this section is that the orlentation 
matrix must be geocentric whether it is referenced to the vernal 
equinox or the Greenwich Meridian. If the orientation matrix is 
geocentric referred to the vernal equinox the geocentric coordinates 
of the camera station are also referred to the vernal equinox. Only 
the end geocentric coordinates of the earth surface point need be 
rotated to be referenced to Greenwich or in the case of a star-referenved 


determination, 
Xx - Xx, cos GST - Ys sin GST 


Y o-= x, sin GST + Ys cos GST 


Tf the orientation matrix is geocentric referenced to 


Greenwich then either GST - 2 replaces 2 in the reduction of 
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Figure 5a GEOCENTRIC ORIENTATION REFERENCED TO GREENWICH 
MERIDIAN OR THE VERNAL EQUINOX 
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Figure 5b ORIENTATION REFERENCED TO LOCAL VERTICAL 
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camera station geocentric direction cosines or the vernal equinox 


geocentric coordinates of the camera station are rotated to Greenwich: 


x = X cos GST - Y sin GST 
fe) fo) fe) 
g s s 

Y = X sin GST + Y cos GST 
re) fe) fe) 
g s s 

Z = Z 
ro) fe) 
g s 


_ Since the most universal form of the orientation matrix 
is star-referenced, let it be assumed that the geocentric coordinates 
of the camera will always be geocentric referred to the vernal equinox. 
Let it be assumed, further, if the orientation matrix is not star- 
referenced that it will be rotated to a star reference whether the 
original matrix refers to the local vertical or Greenwich. Consider 
first a rotation from the local vertical to a star reference. Given 
the orientation matrix referred to the local vertical and the sidereal 


Eulerian angles of the local vertical rotate the matrix to a sidereal 


reference, 
x y Zz 
East (X) cos Ox cos ay cos az 
North (Y) cos Bx cos By cos Bz 


Zenith (Z) cos YX cos Yy cos yz 


x Y Z 
s s s 
East (X) sin 2 cos 2 sin ¢ cos &% cos 
North (Y) cos 2 sin 2 sin ¢ sin 2 cos 6 
Zenith (Z) 0 cos ¢ sin 


where 2 is the right ascension of the local meridian and ¢$ is the 


geodetic latitude defined on the local vertical, 
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Cos OX. = COS OX 
cos dy, = cos aly 
cos AZ. = cos Zz 
cos Bx. = cos ax 
cos BY. = cos ay 
cos Bz. = cos az 
cos YX, = cos 02 
cos Yy, = cos dy 
cos YZ. = cos OZ 


Suppose it is desired to rotate from Greenwich to sidereal. 


x 
Xx cos Ax 
g& 
Y cos Bx 
g B 
Z cos Yx 
g& i & 
xX 
s 
Xx cos GST 
8 
Y - sin GST 
& 
Z 0 
g 


sin 
sin 
sin 
cos 
cos 
cos 
cos 
cos 


cos 


?) 
A?) 
Q 
2 


~~ VD D ®D 


+ cos Bx cos Q 
+ cos By cos 2 
+ cos Bz cos 2 
sin > + cos 8z sin 2 
sin 6 + cos Bz sin 2 
sin 6 + cos Bz sin 2 
cos $¢ + cos £x sin 2 
cos $ + cos By sin Q 


cos ¢ + cos Bz sin Q 


y Zz 
cos 0 cos 0 
Ys “8 
cos B cos Bz 
Ve 8 
cos os 
YV5 c YZ, 
Y Z 
s Ss 
sin GST 0 
cos GST 0 
0 1 


cos OX = cos a cos GST — cos oe sin GST 


cos Oy = cos ce cos GST - cos By, sin GST 
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sin ¢ + 


sin ¢ +c 


sin 6 + 
cos ¢ + 
cos ¢ + 
cos @ + 


cos 
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cos 
cos 
cos o> 
sin 
sin © 


sin 6 
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cos az. = cos Oz cos GST - cos Bz, sin GST 
cos Bx, = cos Ox, sin GST + cos er cos GST 
cos By_ = cos ay, sin GST + cos By, cos GST 
cos Bz. = cos an sin GST + cos Bay cos GST 


cos YX_ = cos yx 
cos Yy_ = cos Vp 
cos YZ, = cos 25 


Clearly if a matrix can be rotated in one direction it can 
be rotated in the opposite or continue in the same direction (360°-GST) 


from sidereal to Greenwich or to the local vertical. 


Henceforth, it will be assumed that the ortentation matrix 
and the coordinates of the camera station are geocentric siderial referenced. 
Then given the camera coordinates of an image of an earth surface point 
the geocentric sidereal reference direction cosines of the line are 


computed: 


cos ox. “i + cos ay Yi + cos az £ = cos a 
eo ae ae? 
i i i 
x Mis 
cos Bx. i + cos BY, i + cos Bz. f = cos B 
ty i ty 
x, y 
cos YX, i + cos eee i + cos YZ, £f = cos Y 
1; 7 i 
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lL, 
2 = 2 2 2 2 
where bs (x; aoe a + £°) 

cos 4, cos 8, cos y of the line, Xo Yo? Zo of the camera 
station, a, b of the earth and h the elevation of the point are the given 
data in the determination of the geocentric coordinates of the earth 


surface point. 


4, THREE METHODS OF DETERMINING THE GEOCENTRIC COORDINATES OF 
AN EARTH SURFACE POINT IN WHICH h fe) 


4.1 Expanded Ellipsoid of Revolution Method 


The geometry of the Expanded Ellipsoid Method is illustrated 


in Figure 6. The equatorial and polar radius of the earth are expanded 
to 


ath 


and 


b' = bith 


producing an enlarged ellipsoid of revolution that places the elevated 
point cofmcident with an imaginary datum. This removes the largest 


error arising from relief displacement in the quadratic determination 
M. 


In the formation of the quadratic 
' 


e2 = (a+ h)* - (b +h)? 
(a + h)- 


whence 


(cos? a + cos* 8) (1 - a?) + cos2y] 
q 


is ig Lee 
2 [x cos a + Xo cos: 6) (1 - e+) + zy cos y] 


w 
i] 
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h COS E (i-e7)% 
(i-e? COS*E) % 


RP=acCoOsE +r. 


A SINE 


=b SINE + —~—~—_— 
é ? (i-e* cos? Ee) 2 


Figure 6 ELLIPSOID a AND b EXPANDED BY h 


26 
Approved For Release 2003/06/20 : CIA-RDP80T00703A000800010001-3 


Approved For Release 2003/06/20 : CIA-RDP80T00703A000800010001-3 


' 1 ' 
= 2 2 un pele 2 _ 2f71 — 2 
c= (xo + YO) Q e“) + Ze a“(1 es) 


and 


Ke 
i] 
N wW 
| l 
Ss & 
QoQ. 6.8 
o Oo 
n if) 
<~ DBD 


Having employed a' and b" to determine the coordinates of 
a point on an expanded ellipsoid of revolution the accepted values of a 
and b are employed to determine latitude (6). Longitude (\) is 


independent of elevation and may be calculated directly: 


tan A =¥ 
x 
sinr> = 
R 
cos A = X 
R 


1 
where R = (X2 + yY2)% 


@ is determined directly or indirectly by an iterative 
process. In these reductions, iteration is employed as a computational 
convenience to avoid the less efficient explicit solution of a fourth 


degree polynomial. 
Parametric angle is E is employed in the indirect method. 


tan E' =Z a 
R b 
27 
Approved For Release 2003/06/20 : CIA-RDP80T00703A000800010001-3 . 


Approved For Release 2003/06/20 : CIA-RDP80T00703A000800010001-3 


where 
R= (xX* + ¥2)% 
Then 
A=e*— [R sec E' - Zb cse E'] 
a aa 
AE = A 7 
a aa 
On the nth iteration 
H= Et + Jaz: 
Where n is the iteration for which 
A=0 
With E determined 
tan d= tan Ea = tanE , 
& GU ee)? 
In the direct method 
tan ¢'= Z at 
R be 
A = -[R sec ¢'(1 - e*%) - Z ese 6' +h e?] 
Ad = A 
Risec 6'(1- e2) tan 6’ + Zese o’ oT 
and on the nth iteration 
d= o'+ ) Ad 
4.1.2 Numerical Example 
Given 
xX, = 8,314,238.0 
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rd 
i 


= 13,305,562.0 


N 
i] 


14,584,184.0 


cos a= .87122760 


cos B= .46326245 


cos y=-. 16233007 


The camera station has an altitude of 100 miles at 


43° 5° 45" 5 


oa 
i] 


58° 0" 0" ,0 


Pg 
n 


The true value of M, ¢, A, and h to be determined are 


as follows: 
M = 1,047,212.0 ft. 
6 = 45° 5' 46",3 


A = 60° 


a 
Il 


10,000.0 ft. 


Initially M is determined 


"= a+t+h = 20,935,689.0 
b' = b +h = 20,865,539.0 
29 
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=a? = b2 = 00669025 


1 - e* = .99330975 
1 
A = {(eos* o + cos2g) (1 - e%) + cos2y] 
= » 99348672 
1 
B = -2[ (cos OX +- cos BY .) (L - e%) + cos yZ 01 
= -21,900,818.0 
1 ’ ' 
7 2 25° V7 nae 2 xo oe eae 
C (X, + YX ) a es) + Zo at (1 e+) 
= 21,845,337 x 10 1? 
M = -B +\/B2 ~ 4ac 
i ) ae 
M = 1,047,215.0 which is only 3 ft, in error. 
whence 
Error 
X = xX, -Mcos a = 7,401,878.0 0.0 
Y= YO - Mcos 8 = 12,820,428.0 0.0 
Z= 2 - Mcos y = 14,754,178.0 0.0 


Thus, there is 


point from the 


no error in the geocentric coordinates of the unknown 


approximation of replacing a and b with a' and b'. 


Next the latitude, longitude are determined and the elevation is 


checked. 
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By the indirect method 


R= (x2 + ¥2)% = 14,803,755.0 


tan E' =Za 
Rb 
= (.99665105) (1.0033,6361) 
= 1.00000340 
BY =. 45° 0° 'O",.35 
sin E' = .70710798 
cos E' = ,70710558 


-R/a sec E' =1.00047872 


esc E' = .99378209 


e2 =-,00669663 


e* = .00669342 


A ==-.-,00000321 


R/a sec E' tan E' = 1,00048212 


Zb ese E' cot E'= "| ,99377871 
aa 1.99426083 


AE =A = -,00000161 


AEt= — 0",33 
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= E' + AE 


teal 
l 


45° 0' 0",35 = 0".33 


45° O' 0",02 
The true value of E is 45° 0’ O".0. Thus, the true value 
is reached to the significant decimal place on the first iteration 


In any case 


tan $= tan Ea ='tan E 


and 
cos A= X = .50000003 
R 
sin A= Y = .86602541 
R 
Error 
o> = 45° 5' 46.334 0" .02 
4 = 60° o"'.0 
Xy = 7,398,3480 0.0 £t 
‘yy =12,814,314.0 -1.0 ft 
Zu =14,747,094.0 -1.0 ft 


N denotes the geocentric coordinates of the point at sea level, 


By the direct method 


32 
Approved For Release 2003/06/20 : CIA-RDP80T00703A000800010001-3 


Approved For Release 2003/06/20 : CIA-RDP80T00703A000800010001-3 


tan ¢' = 


ZAIN 
foal TS) 
NN 


1.00336701 


= 45° 5' 46",.66 (error 0",31) 


oe 
l 


» 70829419 


sin o' 


- 70591738 


lt 


cos o! 


I 
oO 

N 
ed 

tt 


- R sec o' (1 20,830,578,2 
- he? = 6693422 
- [R sec 6' (1 - e*) + he ] = -20,830,645.13442 


+ Z csc $= '20,830,578.8 
= -66.33422 


R sec o' (1 - e”) tan o'= 20,900,782.12 
+ Z esc $' cot $'= 20,760,677.39 
) = 41,661,459,.51 


Ad = = ,00000159 
Ad = -0",33 
6 = o' + Adm 
= 45° 5' 46".33 (error 0,02) 


Thus either method direct or indirect converges with equal efficiency. 
The longitude reduction is the same in either case. A check on the 


reduction is obtained with 


= 
i) 


[Ong Se PE ee ere 


10,001.2 which is an error of 1.2 feet. 


Under the geometric conditions of the model the expanded 
ellipsoid is a simple method yielding accurate results. The relevant 


geometry of the model is h = 10,000 feet, H, = 528,000.0 and ¢, is 61°. 
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GG is the angle enclosed at the camera station between the camera radius 
vector and the line in space. Obviously as on + 0 errors arising from 
relief displacement approach zero. Conversely as ue > 90° the relief 


displacement approaches m- 


4.2 Method Based Zero Change in M, X, ¥ Z and > 


The second method is illustrated in Figure 7. co is the 


angle enclosed between M and P, 


Initially M is co ed as if h were zero using a and b in 
the quadratic in M = B+ \B - 4AC 


OK 


where 


[(cos20 + cos?) (1 - e2) + cosy] 
-2 [ (cos OX + cos BY) (1 - e*) + cos yz! 
[(x,? + yo) G -e*) + Z° - a(l - e”)] 


tb 
It 


Having M to the earth surface p. is the radius vector of this part 


approximated: 


Ul 
Hae 
i 
ik 
Q 
ie) 
a 
g 


tan o' = vA 
er 7 a 
(x2 + Y2)% (1 - e*) 
AR = h sin 9' 
AR = h cos $' 
Zia Z +:AZ 
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Figure 7 METHOD BASED ON ZERO CHANGE IN M, X, Y, Z, AND @ 


35 
Approved For Release 2003/06/20 : CIA-RDP80T00703A000800010001-3 


Approved For Release 2003/06/20 : CIA-RDP80T00703A000800010001-3 


R' =R+ AR 

T t 
ep, = (2% + R4) 
sin ¢, = Bo Sone 50 


Mj = Py cos t = 6, cos v, 


Then approximate values of the geocentric coordinates X, Y, Z and 


latitude (>) are computed: 


Xl x ~ Mi cos a 


a = Y, ~My cos B 


vA =Z,-M, cos y 


Al 
tan $5 = Z, 
R,( - e*) 
or 
tan EY = 41a = al 


e% (81 sec wy - 2 bese BI) 

AE = a a a 
R sec EN tan EY + oS b= cse EY cat EY 
a a 


fon YAR 
Then as in the first iteration 


tan E, 


tan >. = . 
(1 = e*)% 
Bi me 60s Eth cos od. 
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On the 2nd iteration 


X, = Xo - M, cosa = 7,401,967.0 
XY - M, cos 8 = 12,820,475.0 
14,754,161.0 


Z, = Zo - M, cos y 


R, = (X;) + Y¥, ) = 14,803,840.5 


tan E] = Zo 


= 99999647 

Ei = 44° 59" 59.64 
sin E}] = .70710674 
cos E] = ..70710801 


e* —- ["1 sec E] - 4G eh esc EJ] =A = .00000356 
a a 


sec E] tan E} + Ay b csc EI cot El = d= 1.99426494 
a a a 


AE = A = 0,.00000178 
AE" = 0.37 
= ' W 
E, = EY + AE 


= 45° OF O",01 


sin EF, = .70710635 
cos FE. -:= .70710675 


tan E = 1,00336370 
(1 - e2)k-=45° 5' 46'', 32 


ct 

99 

Hs 3 

Ce 
= 

il ll 


14,803,755.0 
14,754,177.0 


R’ = a cos E, +h cos $4 
b sin BE +h sin e, 


N 
i] 


PP Ge f ie: = 20,900,644.0 @. py US ae. 
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= 20,925,689.0 a(1 - e*) = 20,785,625.0 


© 
t 


e2 = ,00669342 


(1 - e* sin*¢,) = .99663737 
(1 - e? sin’o,)% = .99831727 
(1 - e* sin*o,)3/2 = .99496030 
m = 20,900,909.0 
n = 20,970,961.0 
sin >, cos Aj = 0.35394455 cos $, Sin Ay = 0.61117156 
sin >, sin 1 = 0.61408666 cos ¢, cos Ay = 0.35226435 
cos $y = 0.70542246 
Condition Equations 
cos a*AM A+Ao + BeAr = dX 
0.87122760 -7,397,763.0 -12,816,855.0 = -3523.0 
cos §*AM Ao *Ao Bo *AA = dY 
0.46326245 -12,834,969.0 1530/4322 60 = -6112.0 
cos y*AM A3*Ao B3*AA = dZ 
-0,.16233007 14,743,971.0 0.0 = -7088.0 


3 

Let dM = x-10* 
Ad = Ye10 

_4 
Rio ezeLO 
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Zs = b sin E) +h sin o. 


= (R12 12)1; 
Po= (RUS + Z1°)%s 
sin a Po one Fo 


2 
M, = DG cos So - Py cos So 


When two successive values of M, X, Y, Z and » are obtained without 
signficant numerical change final values ¢ and 4 are computed with final 


values of X, Y, Z which are computed with the final value of M. 


4.2.1 Numerical Example 


Using the data from the previous model 


20,925,689.0 
b = 20,855,539.0 
e = 006693422 
(1-e2)= .99330658 
a (l-e”)= 20,785,625.0 


Xo 8,314, 238.0 
¥ = 13,305,562.0 
Zs 14,584,184.0 
cos a = 87122760 
cos B = 46326245 
cos y =-.16233007 
A = E(cos? a + cos28) (1 - e%) + cosy] 
= .99348364 
B = -2[ (cos a x, + cos B XY) (1 ~ e*) + cos y Zo] 
= -21,900,734.0 
c= [x,° + Oe Ci. =84): + z° - a*(1 - e7)] 


22.261,74961 x 107% 
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M = 1,067,893.95 
M'= 20,973,138.51 


The shorter value is the desired value. 


X = Xx, -Mcos a= 7,383,859.0 
Yr Yo - Mcos B = 12,810,847.0 
Z= Zo - M cos y = 14,757,535.0 

R= (X2 + Y2)%= 14,786,452.4 


tan o' = Z 
R(1-e2) 


1.00476968 
o" = 45° 8' 10.74 


AZ = h sin $' 7087.87 
AR = h cos 9' 7054.22 
R+AR=R = 14,793,507.6 

Z+ AZ = Zz = 14,764,623.0 

0 = (R'2 + 2Z'2)% = 20,900,693.0 


= 2 2 2yL, 
Py (XK, + YX + Z \45 
= 21,421,086.73 
cos GG = cos a ae + cos 8B — + cos y iA 
Py Py Po 
= .51538511 


c = 58° 58' 36",43 
sin ¢ = .85695868 
= sin Z 
Py 1 
p 
= .8782,9558 
co = 61° 26' 14",78 
cost, = - 47811806 
P, cos t= 11,040,109.0 


p, cos t= 9,992,999.0 


HN, 1,047,110.0 
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Pern t> 7 Po sin bo 
p 
= .87829764 
ie 61° 26' 15.67 te a 0.89 
cos bo = ,47811427 
= 11.040,109.0 


cos 
Po SO 


5 cos Co a 9,992,896.0 


A 


My = 1,047,213.0 M 
Xy = x, - M, cos a = 7,401,877 ey X, = 90 
Yo = Xo - My cos 8 = 12,820,428 be Y, = 47 
Zo = Zo - M, cos y = 14,754,178 Z,- Zo =~17 
R, = (X,2 + Y,°)% 
2 2 2 
= 14,803,755.0 R, - R= 95.0 
tan E) at SD 
Ry (1 = e*)iy 
= 1,00000341 
EY, a 45° 07 0".35 
cos E, = ,70710558 
sin ES = .70710798 
e? - [82sec BL - “2 = b ese Bh] =A = -.00000320 
a a a 


AL 
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sec BE} tan BE} + “3 b csc E} cot ES =) = 1.99426082 
a a 


a 
Bo 


AE, - A = .00000160 


45° O' 0".02 


Gattis tot che Ey = 1,00336371 
Ean 6 5 


R $5 = 45° 5' 46", 36 by -o, = 07.04 


Ry = a cos E, + h cos 4 = 14,803,754 Roe Ry =. 0 


Z5 = b sin E, +h sin 5 = 14,754,177 12.3 Ly = 1.0 


Since the change in 94, R, and Z is negligible, the change in po, M, and 
S must also be negligible. 


Therefore, the final values are 


D5.™ 45° 5" 46",36 


cos Ay = 2 = .499,999 98 
R 
2 
sin A, = ty .86602544 
R 
2 
Ao = 60 


Latitude converges within 0".04 and the longitude to within 
o".0. 
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4.3 Method Based on Iteration of M, ¢, and A 


In this method M is computed as if h were zero followed by 
adjusting the geocentric coordinates for elevation of the point. The 
geocentric coordinates corrected for h combined with M cos a, Mcos &B, 
M cos y yield Xie Yi» and Zs whose errors are a function of the errors 


in M, ¢, and A, 


Initially > and y are appreximated. 


M= -B +” B*  — 4AC 
2A 


tan $9 = Z 
R (1 - e*)% 
cos A’ = x 
R 
sin A = Y 
R 


Introducing h the geocentric coordinates are adjusted 


a +h cos ¢ COs } 


ni ee 
vl - e2 sin? ¢ 


1 


as 
i] 


KR 
il 
ay) 


+hl cos >? sin 


N 
I 
is) 
oN 
be 

I 
o 
ww 
+ 
ess 
ie) 
fs 
+) 
<a 
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Then 


dX =X - (Mcogga + X): 
1 ° 1 

dY =Y - (Mcosg+ Y) 
1 oO 1 

dZ = Z - (M cos y + 4 


The above residuals are functions of erros in M, ¢, and A. 


Differentiating M, 6, and >} with respect to dx.» dy: and dz. 


dX = - a(l - e7) +h sin ¢ cos A Ad 
(1 - e- sin2 >)3/2 
- a + hi cos ¢ sin A AA 
(1 - e@ sin*® 6)% 
dy = -| a(1 - e”) +h| sin ¢ sin A Ad 
(i — e” sin7¢)3/2 
+ a + h} .cos @ cos A Ad 
(i - e@ sin* do) 
dZ = a(l - e2) + h| cos Ad 
(i - e? sin” $6)3/2 
or 
cos a*AM + A 7 + ae AX = dX 
cos BeAM + Abb + B van = dY 
cos yeAM + A 49 + “BAR = dz 
where 
A =| a(l - e2) + h| sin 6 cos A 
1 (1 - e2 sin2 9)3/2 
B = - a + hj] cos ¢ sin 
Cr - e” sin” 4)% 
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A =-| a (1 - e2) +h sin ¢ sin i 
2 (1 = e? sin? ¢)3/2 

B = a +h cos ¢ cos i 
2 (1 - e2 sin” o )45 

A = a(l - e2) +h cos ¢ 
3 (1 - e@ sin? $)3/2 

B = O 


Solving the above 3 x 3 for AM, Ad, AA 


2 1 1 
bo =o, + 6 
Ayn =A + AA 
200°) 1 
Whence 
Xx = a +h] cos 6 cos x 
2 La — e2 sin7$5 hs 2 
Y, = a + hj} cos ¢ sin 
2 [ae Sint s : : 
25 =| a(l - e7) + hj sin bo 
(1 - e2 sin )3s 
Then again 


dx, = x, - (M, cos a + Xo) 


dy, = XY - (My cos B + Y4) 


dz, 


25 - (My cos y + Zz.) 


again the 3 x 3 is formed and solved for dM. Abo» and AYo« This process is 


repeated until 
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dX = 0 

n 
dy = 0 

n 
dZ = 0 

n 

and 

M = M ~ 1 + dM a 
= -A 
vs De cangt ra OS 
A om Ar 

n a - I if n-l 


are the final values where n is the iteration in which the residuais 


vanish. 
4.3.1 Numerical Example 
Given Xo» Yo Z, and cos a, cos 8, cos y 


M is determined as if h were zero 


M= -B + ¥B~ — 4ac 


2A 
= 1,067,894.0 
= Xo - Mcos a= 7,383,859.0 


XY, -Mcos 8= 12,810,847.0 
ZS -Mcos y= 14,757,535.0 
(x? + y7)% = 14,786,452.0 


Ht 


BN « RK 
tl 


Approximate latitude and longitude 


tan > = Z = 1.00476968 
1 "RT - 62) 
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¢ = 45° 8' 10".74 


sin Ps, = ,70878711 
cos e = 70542246 
i= 1.73497990 
Xx 

The coordinates of the surface point are adjusted for h = 10,000 
ft. 


X, =X +h cos $] cos Ay = 7,387,382.0 


It 


Y, = Y¥+h cos ¢, sin hj 12,816,959.0 


Z 


Z+h sin $4 


14,764,623.0 


dX =X —- (M cos a+ Xj) = -3523.0 


CY cS al. (M cos 8 + Yj) = 6112.0 
dZ = 24 - (M cos y + Z)) = -7088,.0 
Let 

m= a(l - e”) +h 


cee 2 a A a aS Io 
(L — e* sin“¢z)3/2 


a 
Gaeta ee 
(i - e* sin* $4)4s 


A, = -m sin $ cos i BB} = -n cos ¢ sin 
Ao = -m sin > sin x Bo = ncos ¢ cos 
A3 = m cos ¢ B3 = 9 


sin A, = .8663,9084 
cos Ay = -4993,6650 
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Whence 
xX Y Z = r 
871.22760 -739.77630 -1,281.6855 = -3523.0 
463.26245 -1,283.49690 738.7322 = -6112.0 
~162.33007 1,474.39710 0.0 = -7088,0 


Eliminating Z from the first two equations by successive division 


and subtraction: 


(1) 0.67975147 -0.57719019 1.0 = -2.74872424 © 
(2) 0.62710472 -1.73743191 -1.0 = -8.27363420 
(1)&(2) 1.30685619 -2,.31462210 0.0 =-11.02235840 
(3) -162.33007 1,474.3971 0.0 =-7088.0 


Eliminating X from the resultant of (1), (2), and (3): 


(4) 1.0 =1.7711375/ =-8, 43425506 
(5) -1.0 9 .08271092 =-43,66412211 
J sdb o7 330 =-52.09837717 


Y = -7.125,467,34 


ho = Y¥ x 104 = ~,000712550 
Ag'= Ad 
sin 1" 
sin 1" = .00000484813672 
Ag” = -146",97 
= -2' 26,97 


Substituting in (4) and (5) to obtain X in the back solution: 


- 8.43425506 -43,66412211 
-12.62018291 64.71856002 
-21.05443797 +21.05443791 
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1000X 
=21,054.44 ft. 


dM 


Substituting KX and Y in (1) and (2) to obtain Z in the back solution: 


- 2.74872424 - 8,27363420 
- 4.11274985 -12.38001433 
- 6.86147409 -20.65364853 
14.31178514 13.20333741 
2 = 7.45031105 -= 7.45031112 
ee tot 
= ~,.0007450311 
Ak" = Ad/sin 1" 
= -153",67 
m =2' 33.67 


These corrections are applied to the first approximations and 


the condition equations are reformed: 


My 1,067,894 64 AS? OF 9074 AG 60° 2" 30".85 
AM - 21,054 - 2’ 26",97 _ 2' 33", 67 
M, 1,046,840 b5 45% BY OES IE Ks 59° 591 57".18 


In the 2nd iteration 


a +h 


dX, = X - [M, cos a + 
fe) 2 oe 
V1 - e% sin” $5 


cos b>» cos hol 


a +h 


a el 


, T - e% sin $5 


dY 


Xy - IM, cos B + cos b» sin ol 
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cos 


sin 
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= = =) 42 
Zo IM, cos y + a(l e+) + h| sin bo] 
VT - e% sin7$5 
$2 A2 
~ 70592730 ~50001184 
~ 70828429 .86601857 


cos by» cos Ay 0.35297201 


0,61134615 


cos > sin ho 


- 99664213 


ke 
1 
oO 
N 
on 
Hs 
=] 
N 
aa 
Ne 
Wt 


-99831965 


YL — e- sin? o5 


a +h = 20,970,911.0 
a(1l - a2) +h = 20,830,611.0 


xX - [M, cos a + a +h 
fe) 2 ne 
YVI — e- sin? bo 


= 8,314,238 - [912,036.0 + 7,402,145] = 57.0 


cos b>» cos hol = dX 


a +h 


1 - [M, cos 6 + 
VI = e* sin” o 


cos bo sin ol = dV 


13,305.562 - [484,962 + 12,820,486] = 114.0 


a(l - e7) +h 
VI - e* sin 5 


Z - IM, cos A + sin by] = dZ 


14,584,184 = [-169,934 + 14,753,995] = 123.0 
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ii = aU = e2) + h = 20,900,759.0 
(1 - e” sin? ) 3/2 


n= a +h = 20,970,911.0 
(1 - e” sin” bo) 


sin $5 cos \ 0.3541,5053 


2 


sin >» sin i 0.6133, 8735 


2 


Condition Equations 


X Y Z r 
(1) 871,22760 -740, 201488 ~1282.048510 = 57,0 
(2) 463.26245 -1282.026118 740.214461 = 114.0 
(3) -162.33007 1480.367925 0.0 = 123.0 


Solving by successive division and subtraction in the order 


of previous reduction. 


(1) 0.67955897  — .57735838 - 1.0 = 0.04446009 
(2) 0.62584896 _-!1.73196578 1.0 = 0.15400942 
(1) & (2) 1.30540793 = 2,30932416 0.0 = 0.19846951 
(3) -162. 233007 1480.367925 0.0 = 123.0 

(4) 1.0 - 1.769044071 = 0,.15203639 
(5) -1.0 9.12494906 = 0.75816877 
7.35590499 = 0.91020516 

Y= 0.1237380256 

Ao= ¥+10-" = 0,0000123738 

Ad" = Ad/sin 1" = + 02",55 


Back solution for X 


0.15203639 0.75816877 

0.21889802 ~1.12910310 

X 0.37093441 -0.37093441 
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aM = x-10%? = 370.9 ft. 


Back Solution for Z 


04446009 0.15400942 
07144119 0.21431003 
@.11590128 0.36831945 
—- .25207181 - .23214891 
- .13617053 .13617054 
= ze10-+ = ,000013617054 
AA" = Ad/sin 1" = 02".81 
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Summary of Iterations 


M ft. o) A 

lst Approximation 1,067,894 45° 8' 10°74 60° 2' 30°85 
lst Correction - 21,054 - 2' 26797 - 2’ 33°67 
2nd Approximation 1,046, 840 45° 5! 43077 59° 59' 57718 
2nd Correction 371 + 255 + 281 
3rd Approximation 1,047,211 45° 5' 46732 59° 59' 59.99 
True Values 1,047,212 45° 5° 46V31 60° 

Error ~ 1.0 ft. + O'O1 - ov01 


The advantage of the third method is the fact that it converges to the true 


value without the tedious iteration for 


er = B sec ut - 22 ese wt 
ceca ' ' Zb t 1 
— sec E' tan E' +—-— csc E' cot E 
a aa 
or 
- [h,* + R sec $'(l-e*) - Z csc $'] 
Ad = 


~ R sec o'(1-e*) tan 6' + Z ese $' cot ¢' 
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Se NUMERICAL EVALUATION 


The numerical examples indicate the first method is the most 
expeditious and the third method is the most rigorous and general. Con- 
ceivably the expanded ellipsoid method might yield inaccurate position 
coordinates for abnormally high altitudes. The expanded ellipsoid method 
should be tested through the range of altitudes found in the earth surface 
features. A range of h values from 1000 ft. to 30,000 ft. in increasing 
increments of 1000 ft. would cover this geometric aspect of the test. Earth 
surface detail has little or no interpretation value as specific points for 
camera stations higher than 500 miles. Therefore testing the first and 
third methods on earth surface points ranging 1000 ft. to 30,000 ft. at 
camera stations ranging in altitude from 100 miles to 500 miles would provide 
numerical data on all possible geometries that the methods might be applied 
to in practice. A range of 30,000 ft. for h in 1000 ft. increments combined 
with a range of 500 miles for H, in 100 mile increments requires the reduction 
of each of the two methods to 5 x 30 models or 300 reductions. Geocentric 
coordinates of the camera station and the earth surface point and the 
direction cosines of the corresponding lines in space must be generated 
for each model. The latitude ($), longitude (1) and parametric angle (FE) are 
fixed for both the camera station and earth surface point. The geocentric 
coordinates of the camera station (X» Yo? Z,) vary with H; the geocentric 
coordinates of the surface point (X, Y, X)° vary with h; and the direction 


cosines of the corresponding lines vary with either a change in h or H. 


Assume the following constant data. 


Earth Surface Point Camera Station 
E = 45° N ES = 43° N 
A = 60° W do = 58° W 
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tan > = 1,00336361 


The constant geocentric coordinates ar 


N 
i} 


. 70710678 
. 70710678 


cos 


sin 


a 


N 
u 


The geocentric coordinates of the earth surface point 


> = 


45° 


a cos E cos 4 
a cos E sin A 
b sin E 
A 
~ 500000 


. 86602540 


7,398, 348.0 


= 12,814,315.0 


14, 747,095.0 


5° 46°3 


20,925,689 feet 
20,855,539 feet 


«73135370 
-68199836 


station with any value of h and H are simply: 


Ps 
W 


K 
fi 


N 
i] 


The coefficients of 
values of h just as 


sin 5 


Xi, +h cos > cos A x 
ol 

Ye +h cos > sin A tis 

Za +h sin > 251 


h, cos $ cos A, 


the coefficient of H, cos % cos Ao? cos b, 


cos > sin A and sin @ are constant for all 


are constant for all values of H. 


53 


e those 


il 


0.93565171 
= 43° 5° 4505 


at the datum surface. 


a cos E_ cos A 
re) fo) 
a cos E_ sin A 
re) ro) 


b sin E 
oO 


r 
fe) 
.52991926 


. 84804810 


8, 109,927.0 
12,978,596.0 


14,223,443.0 


and camera 


Xk 


+ in » 
Yor H cos o, sin). 


Z ok + H sin $y 


+ H cos >. cos X 
) Oo 


sin and 
fo) 
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$ ¢, 
cos - 70591861 . 73021028 
sin - 70829296 » 68322244 


To provide a numerical example 


h = 10,000 ft. and H = 528,000 ft. 


are reduced. 


h cos 6 cos 4 = 3,530.0 H cos , cos XG = 204,311.0 
h cos > sin A = 6,113.0 H cos by sin 5 = 326,966.0 
h sin = 7,083.0 H sin , = 360,741.0 


whence 


7 x +h cos ¢ cos A= 7,401,878.0 
+h cos > sin »X = 12,820, 428.0 
+h sin %% = 14,754,178.0 


i) 
al 
ond 


These are the quantities to be subsequently determined. 


Pa 
HH 
He 


+ H cos %, cos dG = §8,314,238.0 


fe) ok 
Xo = Yok + H cos $5 sin AG = 13,305,562.0 
=Z, +H sin 6 = 14,584,184.0 
oO ok fo) 


While X,, Y,, 2, are subsequently solved they are employed in the construction 


of the model to determine the exact values of cos a, cos 6, cos Y and M. 


x -X 


cos = 
Bg - ¥, 
fo) 
cos Bf = uM 
ce 
re) 
cos Y = ii 
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M= [(K,- x)? + - 4)? + @- 2,7) 


X - X, = 912,360.0 
0 

Y= ¥, = 485,134.0 
ie) 

Z - Z, = -169,994.0 
QO 


M = 1,047,212.0 ft. 


‘cos 0 = .87122760 


cos 86 - 46326245 


cos Y =-.16233007 


The input to each reduction is Xs Yo Zo? cos a, cos B, cos y and datum 
constants a and b. The values to be determined in each case are M, $ and re 
Any other model is generated by changing h and H as specified. The reason 
for selecting E = 45° is to subject the model to any errors arising from 
the difference in geocentric and geodetic latitude. This difference is 
illustrated in Figure 8. The geocentric latitude is defined by the line 
from the center of the earth to the point and the equator. It decreases 
with h. Geodetic latitude is defined by a perpendicular to an ellipsoid of 
revolution at the point and the equator. It is independent of elevation. 


The difference $ - %. = Ad is zero at the equator and the pole. It can be 


shown that the difference is a maximum at >, = 45°, 


Z 
ten: Oa 
R(1-e*) 
Z 
tan >. = R 
tan 
tan > = c 
(1-e*) 
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re 8 GEOCENTRIC AND GEODETIC LATITUDE 


Figu 
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tan > - tan %. 
tan AO = 7T fan 6 tan o, 


Substituting 


2 
tan bo 


* (de?) 


1 


tan >. - (1-e?)tan %. 


= I 


(l-e*) + tan*d, 


2 
ie tan o. 


(l-e7) + tan’ | 


Be 
e“sin o. cos >. 


a mmmennemeal 


Done 
1 - e“cos >. 


2 
e“sin 20, 


= TE, 


2(1-e7 cos”) 


sin 26 is a maximum when 
c 
o. = 45° 


and when 


o = 45° 


sin o = cos o, 


2 
e ; 
l-e*cos*>, =]- 3 «Sin 20, 


r] 

re 

1 
ne 


eS 
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or tan is max when Ag | = 45" and when $. = 45° 


tan Ado= e 


= _+0066934 


- 003357938 


Ag" tan Ad 


“ain 


u 


692" 
11" 32" 


Thus the choice E = 45° subjected the methods to the maximum 


difference between geocentric and geodetic, 
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or tan Ad is maximum when 


and when 


> = 45° 


tan Ad 


- 0066934 
1.9933066 


.003357938 


tan Ad 


Tt —_ 
aoe sin 1" 


ues Pa 


Thus, the choice E = 45° subjects the methods to the maximum 


difference between geocentric and geodetic. 
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